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Abstract 

Probabilistic solutions of the so called Schrodinger boundary data problem 
provide for a unique Markovian interpolation between any two strictly posi- 
tive probability densities designed to form the input-output statistics data for 
the process taking place in a finite-time interval. The key issue is to select 
the jointly continuous in all variables positive Feynman-Kac kernel, appro- 
priate for the phenomenological (physical) situation. We extend the existing 
formulations of the problem to cases when the kernel is not a fundamen- 
tal solution of a parabolic equation, and prove the existence of a continuous 
Markov interpolation in this case. Next, we analyze the compatibility of this 
stochastic evolution with the original parabolic dynamics, while assumed to 
be governed by the temporally adjoint pair of (parabolic) partial differential 
equations, and prove that the pertinent random motion is a diffusion process. 
In particular, in conjunction with Born's statistical interpretation postulate in 
quantum theory, we consider stochastic processes which are compatible with 
the Schrodinger picture quantum evolution. 



I. Motivation: Schrodinger's interpolation prob- 
lem through Feynman-Kac kernels 

The issue of deriving a microscopic dynamics from the (phenomenologically or nu- 
merically motivated, by approximating the frequency distributions) input-output 
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statistics data was addressed, as the Schrodinger problem of a probabilistic inter- 



polation, in a number of publications |lfl-[10]. We shall consider Markovian propa- 



gation scenarios so remaining within the well established framework, where for any 
two Borel sets A, B C R on which the respective strictly positive boundary densities 
p(x, 0) and p(x, T) are defined, the transition probability m(A, B) from the set A 
to the set B in the time interval T > has a density given in a specific factorized 
form: 

m(x,y) = f(x)k(x,0,y,T)g(y) 



m(A,B) = dx dym(x } y) 
J A Jb 

dym(x, y) = p(x, 0) , J dxm(x, y) = p(y, T) (1) 

Here, f{x),g(y) are the a priori unknown functions, to come out as solutions 
of the integral (Schrodinger) system of equations (1), provided that in addition to 
the density boundary data we have in hands any strictly positive, continuous in 
space variables function k(x,0,y,T). Our notation makes explicit the dependence 
(in general irrelevant) on the time interval endpoints. It anticipates an important 
restriction we shall impose, that k(x, 0, y, T) must be a strongly continuous dynami- 
cal semigroup kernel: it will secure the Markov property of the sought for stochastic 
process. 

It is the major mathematical discovery M that, without the semigroup assump- 
tion but with the prescribed, nonzero boundary data p(x, 0), p(y, T) and with the 
strictly positive continuous function k(y,0,x,T), the Schrodinger system (1) of in- 
tegral equations admits a unique solution in terms of two nonzero, locally integrable 
functions f(x),g(y) of the same sign (positive, everything is up to a multiplicative 
constant). 

If k(y,0,x,T) is a particular, confined to the time interval endpoints, form of 
a concrete semigroup kernel k(y,s,x,t),0 < s < t < T, let it be a fundamental 
solution associated with (5) (whose existence a priori is not granted), then there 
exists [§, §, |, |K], 0] a function p(y, s, x, t): 

p{y, s, x, t) = k(y, s, x, t) ^ - (2) 

where 

6(x,t) = dyk{x,t,y,T)g(y) (3) 



Q*(y,s)= / dxk(x,0,y,s)f(x) 
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which implements a consistent propagation of the density p(x,t) = 6(x,t)8*(x,t) 
between its boundary versions, according to: 

p(x,t) = J p(y,s,x,t)p(y,s)dy (4) 

< s < t < T 

For a given semigroup which is characterized by its generator (Hamiltonian), the 
kernel k(y,s,x,t) and the emerging transition probability density p(y,s,x,t) are 
unique in view of the uniqueness of solutions f(x),g(y) of (1). For Markov processes, 
the knowledge of the transition probability density p(y, s, x, t) for all intermediate 
times < s < t < T suffices for the derivation of all other relevant characteristics. 

In the framework of the Schrodinger problem the choice of the integral kernel 
k(y,0,x,T) is arbitrary, except for the strict positivity and continuity demand. As 
long as there is no "natural" physical motivation for its concrete functional form, 
the problem is abstract and of no direct physical relevance. 

However, in the context of parabolic partial differential equations this " natural" 
choice is automatically settled if the Feynman-Kac formula can be utilized to repre- 
sent solutions. Indeed, in this case an unambigous strictly positive semigroup kernel 
which is a continuous function of its arguments, can be introduced for a broad class 
of (admissible [|ll|) potentials. Time dependent potentials are here included as well 
|12| , Moreover, in Ref. || we have discussed a possible phenomenological sig- 
nificance of the Feynman-Kac potentials, as contrasted to the usual identification of 
Smoluchowski drifts with force fields affecting particles (up to a coefficient) in the 
standard theory of stochastic diffusion processes. 

In the existing probabilistic investigations [|, [TJ], [8], |S| , based on the exploita- 
tion of the Schrodinger problem strategy, it was generally assumed that the kernel 
actually is a fundamental solution of the parabolic equation. It means that the ker- 
nel is a function with continuous derivatives: first order-with respect to time, second 
order-with respect to space variables. Then, the transition probability density de- 
fined by (2) is a fundamental solution of the Fokker-Planck (second Kolmogorov) 
equation in the pair x, t of variables, and as such is at the same time a solution of the 
backward (first Kolmogorov) equation in the pair y, s. This feature was exploited 
in 

There is a number of mathematical subtleties involved in the fundamental solu- 
tion notion, since in this case, the Feynman-Kac kernel must be a solution of the 
parabolic equation itself. In general, Feynman-Kac kernels may have granted the 
existence status, even as continuous functions [11], [13], |15[, but may not be differen- 
tiable, and need not to be solutions of any conceivable partial differential equations. 
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To our knowledge, this complication in the study of Markovian representations 
of the Schrodinger interpolating dynamics (and the quantum Schrodinger picture 
dynamics in particular) has never been addressed in the literature. Moreover, it 
is far from being obvious that this Markovian interpolation actually is a diffusion 
process. 



II. Schrodinger's interpolation problem: general 
derivation of the stochastic evolution 

1. The Schrodinger system of integral equations 

We shall complement our previous analysis [§, |§ by discussing the issue in more 
detail. It turns out the the crucial step lies in a proper choice of the strictly positive 
and continuous function k(y, s,x,t),s < t which, if we want to construct a Markov 
process, has to satisfy the Chapman-Kolmogorov (semigroup composition) equation. 
To proceed generally let us consider a pair of partial differential equations for real 
functions u(x,t) and v(x, t): 

d t u(x, t) = Au(x, t) — c(x, t)u(x, t) (5) 

d t v(x, t) = —Av(x, t) + c(x, t)v(x, t) 

where, we have eliminated all unnecessary dimensional parameters. 

Usually, [T5|, |TT|, c(x,t) is assumed to be a continuous and bounded from below 
function. We shall adopt weaker conditions. Namely, let us decompose c(x, t) into a 
sum of positive and negative terms: c(x,t) = c + (x,t) — c^(x,t) , c± > where (a) 
C-(x,t) is bounded, while (b) c + (x,t) is bounded on compact sets of R x [0, T]. It 
means that c(x, t) is bounded from below and locally bounded from above. Clearly, 
c(x, t) needs not to be a continuous function and then we encounter weak solutions 
of (5) which admit discontinuities. 

With the first (forward) equation (5) we can immediately associate an integral 
kernel of the time-dependent semigroup (the exponential operator should be under- 
stood as the time-ordered expression): 

k{y,s,x,t) = [exp{- f H{T)dr)]{y,x) (6) 

J s 

where H(t) = — A + c(r). It is clear, that for discontinuous c(x, t), no fundamental 
solutions are admitted by (5). 
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By the Feynman-Kac formula, |13], I2f , we get 

k(y,s,x,t) = I exp[- f c{u(T),T)dT}dfi i (l'f ) (uj) (7) 



where dfi^t) { u ) is the conditional Wiener measure over sample paths of the standard 
Brownian motion. 

It is well known that k is strictly positive in case of c(x, t) which is continuous 
and bounded from below; typical proofs are given under an additional assumption 
that c does not depend on time |n|. However, our assumptions about c(x,t) were 
weaker, and to see that nonetheless k is strictly positive we shall follow the idea 
of Theorem 3.3.3 in |TjJ. Namely, the conditional Wiener measure ^fot) can be 



written as follows 

d$2) = IMt ' s)}- l/2 exp[- { ^-^-} du$> (8) 



4(t - s) 

where dv 1 ^ 8 ) is the normalised Wiener measure JTT[ . We can always choose a certain 
number r > to constrain the event (sample path) set 

Vt(r) = [u : X s (u) = y,X t (uj) = x, sup s < T < t \X T (u)\ < r] (9) 

It comprises these sample trajectories which are bounded by r on the time interval 
[s, t}. In the above, X t (u) is the value taken by the random variable X(t) at time t, 
while a concrete u-th path is sampled. By properly tuning r, we can always achieve 

> \ (10) 
which implies that 

k(y, s, x, t) > i[47r(* - s)]- 1/2 exp[-^-^-} exp[-(t - s)C] > (11) 
2 4(£ — s) 

C = sup s < T <t,uen(r) c+(X t (u),t) 

where, by our assumptions, c + is bounded on compact sets. Consequently, the kernel 
k is strictly positive . 

With the Schrodinger boundary data problem on mind, we must settle an issue 
of the continuity of the kernel. To this end, let us invoke a well known procedure 



of rescaling of path integrals [TT|, |16[: by passing from the "unsealed" sample paths 
u(t) over which the conditional Wiener measure integrates, to the "scaled" paths of 
the Brownian bridge, the (y, x) conditioning can be taken away from the measure. 
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Then, instead of sample paths u connecting points y and x in the time interval 
t — s > 0, we consider the appropriately "scaled" paths of the Brownian bridge a 
connecting the point with again, in the (scaled) time 1. It is possible, in view of 



the decomposition [11, 16 



T 



u(t) = ( )y + ( k+Vt-sa( ) (12) 

v ; K t-s t-s 1 y t-s t-s 1 y t-s t-s 1 y ' 

where a stands for the "scaled" Brownian bridge. Then, we can write 

k(y,s,x,t) = [^(t - s ))-y^exp[-^^-) J d/i(a). (13) 

r f* ~ T T — S , T — S. .,, 

exp\- c(- y + x + Vt - s a(- ) ,T)dr\ 

J s t — s t — s t — s 

where dfi(a) = d^o'iu^) is the normalized Wiener measure integrating with respect 
to the "scaled" Brownian bridge paths, which begin and terminate at the origin 
in-between "scaled time" instants: corresponding to r = s and 1 corresponding to 

T = t. 

This representation of k, if combined with the assumption that c(x, t) is a contin- 
uous function, allows to conclude, [TT|| , that the kernel is continuous in all variables. 
However, our previous assumptions were weaker, and it is instructive to know that 
through suitable approximation techniques, Theorem B.7.1 in Ref.|T^[ proves that 
the kernel is jointly continuous in our case as well. 

It is also clear that k(y,s,x,t) satisfies the Chapman-Kolmogorov composition 
rule. So, the first equation (5) can be used to define the Feynman-Kac kernel, 
appropriate for the Schrodinger problem analysis in terms of a Markov stochastic 
process. 

Let us consider an arbitrary (at the moment) pair of strictly positive, but not 
necessarily continuous, boundary densities po(x) and pr(x). By Jamison's principal 
theorem ||] there exists a unique pair of strictly positive, locally (i.e. on compact 
sets) integrable functions f(x) and g(x) solving the Schrodinger system (1), e.g. 
such that po(x) = f(x) J k(x, 0, y, T)g(y)dy and p T (x) = g(x) J k(y,0, x,T)f(y)dy 
with the kernel k(y,s,x,t) given by (7). 

Let us define: 

g(x,t) = J k(x,t,y,T)g(y)dy , f(x,t) = jk(y,0,x,t)f(y)dy (14) 

The above integrals exist at least for almost every x so that there appears the 
problem of the existence of a unique and continuous transition probability density 
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p(y, s,x,t), (2). We shall assume that the function g(y) is bounded at infinity. This 
means that there exists a constant C > and a compact set K C R such that 
g(y) < C for all y G R\K. Then, for alH < T and any sequences h n — > 0, s n — > 0, 
as n — > oo we get (Zim stands for iim^-^): 

/im|5((x + /i n ,t + s n )-5f(a;,t)| < lim\ Jjk(x + h n ,t + s n ,y,T)-k(x,t,y,T)]g(y)dy\ + 
lim\ [k(x + h n ,t + s n ,y,T)-k(x,t,y,T)}g(y)dy\< (15) 

JR\K 

Km S u Ps€K \k( X + fc„, * + S „, », D - *<*, t, », T)| + 

C-lim \ \k(x + h n ,t + s n ,y,T) - k(x,t,y,T)\dy 
Jr\k 

The first term tends to zero because k is jointly continuous and g is locally in- 
tegrate. The second one tends to zero because of the Lebesgue bounded conver- 
gence theorem. Consequently, our assumption suffices to make g(x,t) continuous on 
R x [0, T). Similarly, we can prove that g(x,t) is bounded. 

Now, we can set according to (2), p(y,s,x,t) = k(y,s,x,t)g(x,t)/g(y,s). Then, 
p(y,s,x,t), < s < t < T becomes a transition probability density of a Markov 
stochastic process with a factorized density p(x,t) = f(x,t)g(x,t). Clearly, this 
stochastic process interpolates between the boundary data po an d pr as time con- 
tinuously varies from to T. Notice that (15) implies the continuity of p in the time 
interval [0, T). 

Although p(y, s, x, i) is continuous in all variables, we cannot be sure that the 
interpolating stochastic process has continuous trajectories, and no specific (e.g. 
Fokker- Planck) partial differential equation can be readily associated with this dy- 
namics. Therefore, we must explicitly verify whether the associated process is 
stochastically continuous. If so, we should know whether it is continuous (i.e. ad- 
mits continuous trajectories). Eventually, we should check the vailidity of conditions 
under which the investigated interpolation can be regarded as a diffusion process. 
The subsequent analysis will prove that this ultimate goal results only due to the 
gradual strengthening of conditions imposed on the parabolic system (5). 

2. Stochastic continuity of the process 

Apart from the generality of formulation of the Schrodinger interpolation problem 
which appears to preclude an unambigous identification (diffusion or not) of the 
constructed stochastic process, we can prove in the present fundamental 
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property of a stochastic dynamics called a stochastic continuity of the process. In 
this connection, compare e.g. || |17| and |18|], where this property is linked to 
the uniqueness of the corresponding Markov semigroup generator. The stochastic 
continuity property is a necessary condition for the process to admit continuous 
trajectories. 

The stochastic process is stochastically continuous, if for the probability of the 
occurence of sample paths u, such that the random variable values X t (uj) along the 
trajectory obey \X t {u) — X s (uS)\ > e , s < t, the following limiting behaviour is 
recovered 

lim tis P[u : \X t (oj) - X s (u)\ > e] = (16) 

for every positive e. This demand can be written in a more handy way in terms of 
the transition probability density p(y, s, x, t) and the density p(x, t) of the process: 

/+oo r 
dyp(y,s) p(y,s,x,t)dx] = (17) 

-oo J\x— y\>e 

So, for the transition density to be stochastically continuous, it suffices that 

lim^siQ / p(y, s, x, s + As)dx = (18) 

J\x-y\>e 

for almost every y e R. 

In view of our construction, (2), we have: 



lim^ s io / p(y, s, x, s + As)dx = (19) 

J\x-y\>e 

r r+co 

/ dx k(y, s, x, s + As) / k(x, s + As, z, T) g(z)dz 

J\x— v\>e J — oo 



l\x-y\>e 

1 

g(y, S) J\x-y\>e 

By changing the order of integrations (allowed by positivity of the involved functions) 



we get: 

r X r+oo 

lim Asl0 p(y, s, x, s + As)dx = — rliniAsio dz g(z) (20) 

J\x-y\>e g{y,S) J-oo 

[ / dx k(y, s, x, s + As) k(x, s + As, z, T)] 

J\x-y\>0 

Because the potential is bounded from below, c > — M for some M > 0, we easily 
arrive at the estimates (use the "scaled" Brownian bridge argument) 

k(y, s,x,s + As) < (AnAs)- 1/2 exp[- ^ X ~^ ] exp(MAs) (21) 
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and 



k(x, s + As, z, T) < [4tt(T - s - As)]~ 1/2 exp[ ! ——}exp[M{T 

4(T — s — As) 

Then we get: 

< lim&sio / k(y, s,x, s + As)k(x, s + As, z, T)dx < 

J\x-y\>e 

[4tt(T - s)}- 1/2 exp[M(T - s)}lim Asl0 (AnAs)- 1/2 (23) 

f i r ( x ~ V) 2 i r ( Z ~ X ) 2 i 

l- yl J XeXp[ -^As- ] eXp[ ~ A{T - s - A S ) ] = ° 

So, by the classic Lebesgue bounded (dominated) convergence theorem, the required 
limiting property lim& s ^ J\ x _ y \ >e p(y, s, x, s+As)dx = follows and (16) holds true. 

As mentioned before, the stochastic continuity of the Markov process is a nec- 
essary condition for the process to be continuous in a more pedestrian sense, i. e. 
to admit continuous sample paths. However, it is insufficient. Hence, additional 
requirements are necessary to allow for a standard diffusion process realization of 
solutions of the general Schrodinger problem, (l)-(3). 

In the next section we shall prove that our process can be regarded as continuous, 
by requiring a certain correlation between the kernel k(y,s,x,t) and a function 
g(x,t), (14). 



- s - As)] 
(22) 



3. Continuity of the process 

It is well known that a solution of a parabolic equation cannot tend to zero arbi- 
trarily fast, when \x\ — ► oo, [|i"9"| . Roughly speaking, it cannot fall off faster than a 
fundamental solution (provided it exists). In fact, the solution is known to fall off 
as fast as the fundamental solution, when the initial boundary data coincide with 
the Dirac measure. If a support of the initial data is spread (i.e. not point-wise), 
then the solution falloff is slowlier than this of the fundamental one. 

In our discussion, where g(x,t) is a generalized solution and k(y,s,x,t) is a 
Feynman-Kac kernel which does not need to be a fundamental solution, we expect 
a similar behaviour. Mathematically, our demand will be expressed as follows. Let 
t — s be small and K be a compact subset in R. Because g(x, t) is supported on the 
whole R, so in the decomposition 

g(y,s)= k(y,s,x,t)g(x,t)dx + k(y, s,x,t)g(x,t)dx (24) 

JK JR\K 
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the second term becomes relevant when \y\ — > oo . It amounts to (in the denominator 
there appears g(y,s)): 



k(y, s, x, t)g(x, t)xi<{x)dx 
k(y,s,x,t)g(x,t)dx 



hm \y\^oo r+ oo ,w._ = ( 25 ) 



where xk is an indicator function of the set K, which is equal one for x G K and 
zero otherwise. 

By means of the transition probability density p(y,s,x,t) let us introduce a 
transformation 

r+oo 

(T t J)(y)= p(y } s } x } t)f(x)dx (26) 
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of a function f(x), continuous and vanishing at infinity (we shall use an abbreviation 
/ ^ Coo(R) to express this fact). It is clear that (T*f)(x) is a continuous function. 
For a suitable compact set K we can always guarrantee the property \f(x)\ < e for 
every x G R\K. Then, if we exploit the property j R \ K p(y,s,x,t)dx < 1 if s < t 
and the definition of p in terms of k and g, we arrive at 



\(T!f)(y)\ < / p(y,s,x,t)\f(x,t)\dx + / p(y,s,x,t)\f(x,t)\dx < (27) 
Jk Jr\k 



p(y,s,x,t)dx] / \f(x,t)\dx + sup xeR \ K \f(x,t)\ / p(y,s,x,t)dx < 

K JK JR\K 

\f(x,t)\dx] lKk{y,s,x,t)g{x,t)dx + e 
^ ' J-™k(y,s } x,t)g(x,t)dx 

It implies that for small t — s, lim\ y \^ QO {T t s f){y) = 0, and so T s * forms an inhomo- 
geneous in time semigroup of positive contractions on Coq(R). For arbitrary t and 
s the result follows by the obvious decomposition property T* = T^T^ ■ ■ ■ T* n . In 
the well established terminology, our p(y, s, x, t) is a Coo-Feller transition function 
and leads to a regular Markov process, [17] . Moreover, by the stochastic continuity 
of p(y, s, x, t), T* is strongly continuous. 

As yet, we do not know whether the process itself is continuous i.e. has contin- 
uous sample paths. To this end, it suffices to check whether the so called "Dynkin 
condition", [pOf 



lim tls —^—sup yeK [[ p(y,s,x,t)dx] = (28) 

t — S J \x— j/|>e 

is valid for every e > and every compact set K. We have (remember that g(x,t) 
is strictly positive, continuous and bounded): 

sup yeK [ p(y,s,x,t)dx = sup yeK , -, / k(y,s,x,t)g(x,t)dx < 

J\x-v\>€ Q(V,S) J\x-v\>e 
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sup x g(x,t) 



/ k(y,a,x,t)dx < C k (x — y,t — s)dx (29) 

J\x—y\>e J\x— y\>e 



infye K g(y } a) 
where (compare e.g. the previous estimate (22)) 



sup x g{x,t) 

C = — r exp[M(t - s) (30) 

mfyeK9{y,s) 



and ko(x — y,t — a) is the heat kernel. 
Finally, we arrive at: 

lim tls —^— aup yeK [ [ p(y, s, x, t)dx] < 

t — S J\x-y\>e 

C lim t is / k (z,t — s)dz = (31) 

t — S J\z\>e 

So, the stochastic process we are dealing with, is continuous. Interestingly, "a con- 
tinuous in time parameter stochastic processes which possesses the (strong) Markov 
property and for which the sample paths X(t) are almost always (i.e. with proba- 
bility one) continuous functions of t is called a diffusion process", see e.g. chapter 
15 of 1201. 



4. The interpolating stochastic dynamics: compatibility with 
the temporally adjoint parabolic evolutions 

The formulas (14) determine what is called, |12j| , the generalized solution of a 
parabolic equation: it admits functions which are not necessarily continuous and 
if continuous, then not necessarily differentiate. Before, we have established the 
continuity of the generalized solution g(x,t) under rather mild assumption about 
the behaviour of g(x) at spatial infinity. In fact, the same assumption works for 
f(x,t). But nothing has been said about the differentiability of f{x,t) and g(x,t). 

Consequently, our reasoning seems to be somewhat divorced from the original 
partial differential equations (5), for which we can take for granted that certain 
solutions u(x,t) and v(x,t) exist in the time interval < t < T. For this, we must 
assume that c(x, t) is a continuous function. 

Let us consider the solutions of (5) that are bounded functions of their argu- 
ments. It is instructive to point out that we do not impose any restrictions on the 
growth of c(x, t) when \x\ — > oo, and consequently we do not assume that solu- 
tions of parabolic equations (5) have bounded derivatives. Then, [O, the solution 
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u(x, t) of the forward parabolic equation (5) is known to admit the Feynman-Kac 
representation with the integral kernel (7), (13), where 



u(x, t) = / k(y, s, x, t)u(y, s)dy (32) 

for < s < t < T. At this point let us define 

U(x,t) — v(x,T — t) (33) 

for all t G [0, T] and observe that, as a consequence of the time adjoint equation 
(5) for which v(x, t) is a solution, the newly introduced function U(x,t) solves the 
forward equation (5): 

d t U(x, t) = AU(x, t) - c(x, T - t)U(x, t) (34) 

with a slightly rearranged potential: c(x,t) — > c(x,T — t). By the assumed bound- 
edness of the solution v(x,t) of (5), we arrive at the Feynman-Kac formula 

U(x, t)= f K{y, s, x, t)U(y, s)dy (35) 



with the corresponding kernel K(y, s, x, t) of the (time ordering implicit) operator 
exp[— f 8 H(T — r)dr], where H(T — r) = —A + c(T — r). Let us emphasize that in 
case of the time independent potential, c(x, t) = c(x) for all < t < T, the kernel 
K coincides with k. 

The previous Brownian bridge argument (12), (13) retains its validity, and we 
have: 

K(y, s, x, t) = [47r(t - a )]-Wexp[-&0£]- (36) 

J dfi(a) exp[— J c(- y + -x + y/t — s a(- -) , T — r)dr] 

which, after specializing to the case of s = 0, t = T and accounting for the invariance 
of the Brownian bridge measure with respect to the replacement of sample paths 
u){t) by sample paths uj{T — r), |7|, [27], gives rise to: 

K(y,0,x,T) = (AnTr^expl-^^]. (37) 

J d/i(a)exp[- J o c(^y + (1 - ^)x + VT a(^),a)da] 
where a = T — t. 
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A comparison of (37) with (13) proves that we have derived an identity: 



K(y,0,x,T) = k(x,0,y,T) (38) 
whose immediate consequence is the formula 

U{x, T) = v(x, 0) = ( k{x, 0, y, T)v(y, T)dy (39) 



for the backward propagation of v(y,T) into v(x, 0). 

We shall utilize (39) and (32), under an additional assumption that the previous, 
hitherto arbitrary, probability density data p (x), px(x), actually are determined by 
the initial and terminal values of the solutions u(x,t), v(x,t) of (5), according to: 

Po(x) = u(x, 0)v(x, 0) 

p T (x) = u{x,T)v{x,T) (40) 

Our present aim is to show that with this assumption, we can identify the (still 
abstract) functions f(x,t), g(x,t), (14), with u(x,t) and v(x,t) respectively. By 
(32), (39) there holds: 

Po(x) = u{x, 0) / k{x, 0, y, T)v{y, T)dy (41) 



p T (x) =v(x,T) / k(y,0,x,T)u(y,0)dy 



and, in view of the uniqueness of solution of the Schrodinger system, once the bound- 
ary densities and the continuous strictly positive kernel are specified, we realize that 
the propagation formulas (14) involve solutions of (5) through the respectively initial 
and terminal data: 

/(a?) = u(x, 0) 

g(x)=v(x,T) (42) 

Moreover, (5), (14) imply that f(x, t) = u(x, t) holds true identically for all t G [0, T]. 

What remains to be settled is whether the function g(x, t) can be identified with 
the solution v(x, t) of (5) for all t G [0,T]. 

This property is obvious, when the time independent potential c(x) is inves- 
tigated instead of the more general c(x,t). As well, the identification is with no 
doubt in case when k(y,s,x,t) is a fundamental solution of the parabolic equation 
in variables x, t. In this case, k(y, s, x, t) is a unique solution of the system (5), and 



solves the adjoint equation in variables y, s, [^TJ |22], [23|. Then, because f(x),g(x 



are locally integrable, an immediate consequence is, that f(x,t) and g(x,t) 
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are positive solutions of (5). The identification of them with u(x,t) and v(x,t) 



respectively, follows from the uniqueness of positive solutions, [22 



Let us begin from a minor generalization of (22), and define: 

U a (x,t)=v{x,T + s-t) , te[s,T] (43) 

Clearly, a parabolic equation (34) is satisfied by U s (x, t), if instead of c(x, T — t), the 
potential c(x, T + s — t) is introduced. An immediate propagation formula follows 

U s (x, t) = J K s (y, s, x, t)U s (y, s)dy (44) 

The integral kernel K s differs from the previous K, (36), in the explicit time depen- 
dence of the potential c(x, T — r) — ► c(x, T + s — r). By putting T = t in (44) we 
get: 

v(x,s)= / K s (y,s,x,T)v(y,T)dy (45) 



and by the previous part of our demonstration we know that 

g(x, s) = J k(x, s, y, T)v(y, T)dy (46) 

At this point, it is enough to prove that the identity (cf. (38)) 

K s (y, s, x, T) = k(x, s, y, T) (47) 

takes place for any s; < s < T . 

Let us exploit the Brownian bridge scaling (13) again, so that 

k(x, a , y, T) = IMT - s)]- 1/2 exp[-^^r}- (48) 

/ rp rp ^_ ^_ ^ ^_ ^ 
dfi{a)exp[- / c(- x + y + VT - s a(- ), r)dr} 
Js 1 — S 1 — S 1 — S 

and, analogously 

K.(y,8,x,T) = [A-k(T - s)]-^exp[-^p^r]- (49) 

4(T - s) 

/ dfi(a)exp[- / c(- y + x + VT - s a(- ), T + s- r)dr] 

J Js 1 — S 1 — S 1 — S 

By changing: 

«(^) =► ocil - 1^) = a(^) (50) 
1 — s 1 — s 1 — s 
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and substituting a = T + s — r, where r only is the running variable, we finally 



recover 



K s (y,s,x,T) = [47r(T- S )]-V2 ea ;p[_fe__^!]. ( 5 1) 
f 1 \ r r i a ~ s T — a r— o — s 

/ dfi{a)exp[- / c(- y + - x + VT - s a{- ), a){-da)\ = k{x,s,y,T) 

J JT 1 — S 1 — S 1 — S 

Hence, 

g(x, s) = v(x, s) (52) 

is valid for all time instants < s < T. This implies that p(y, s, x, t) = k(y, s, x, t) ^r'x 
defines a consistent transition probability density of the continuous Markovian in- 
terpolation. 



We have succeeded to prove that: 



(i) If a continuous, strictly positive Feynman-Kac kernel of the forward parabolic 
equation (5) is employed to solve the Schrodinger boundary data problem (1) for an 
arbitrary pair of nonzero probability densities po(x) and Pt{x), then we can construct 
a Markov stochastic process, which is continuous and provides for an interpolation 
between these boundary data in the time interval [0,T]. 



(ii) Given the time adjoint parabolic system (5) with bounded solutions u(x,t), 
v(x,t) in the time interval [0, T\. If the boundary densities are defined according 
to (40), then the Schrodinger problem (l)-(3) provides us with a unique continu- 
ous Markov interpolation, that is compatible with the time evolution of p(x, t) = 
u(x,t)v(x,t), t G [0,T]. 



5. Whence diffusions ? 

Our strategy, of deducing a probabilistic solution of the Schrodinger boundary data 
problem in terms of Markov stochastic processes running in a continuous time, was 
accomplished in a number of steps accompanied by the gradual strengthening of 
restrictions imposed on the Feynman-Kac potential, to yield a continuous process 
(cf. Section II. 3), and eventually to get it compatible with a given a priori parabolic 



evolution (Section II. 4). In a broad sense, [ 20fl , it can be named a diffusion. 

However, this rather broad definition of the diffusion process is significantly 
narrowed in the physical literature: while demanding the continuity of the process, 
the additional restrictions are imposed to guarrantee that the mean and variance 
of the infinitesimal displacements of the process have the standard meaning of the 
drift and diffusion coefficient, respectively, 
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According to the general wisdom, diffusions arise in conjunction with the parabolic 
evolution equations, since then only the conditional averages are believed to make 
sense in the local description of the dynamics. It is not accidental that forward 
parabolic equations (5) are commonly called the generalized diffusion equations. 
Also, the fact that the Feynman-Kac formula involves the integration over sam- 
ple paths of the Wiener process, seems to suggest some diffusive features of the 
Schrodinger interpolation, even if we are unable to establish this fact in a canonical 
manner. 

Clearly, the conditions valid for any e > 0: 

(a) there holds lim t i s jz^ J\ y _ x \ >€ p(y, s,x,t)dx = 0, (notice that (a) is a direct con- 
sequence of the stronger, Dynkin condition, (28)), 

(b) there exists a drift function b(x, s) = lim t i s -^ Ju_ a i< e (y — x)p(x, s,y,t)dy, 

(c) there exists a diffusion function a(x, s) = Um t i s -^ J\ y - X \<e(y ~ x) 2 p(x, s, y, t)dy, 
are conventionally interpreted to define a diffusion process, ]26| . 

To our knowledge, no rigorous demonstration is available in the Schrodinger 
problem context, in case when the involved semigroup kernel is not a fundamental 
solution of the parabolic equation. 

Let us impose a restriction on a lower bound of a solution v (x, t) of the backward 
equation (5). Namely, we assume that there exist constants c\ > 0, c 2 > such that 
v(y,s) > Ciexp(—C2y 2 ) for all s e [0,£],t < T. This property was found to be 



respected by a large class of parabolic equations, [25], and it automatically ensures 



that the condition (25) of Section II. 3 is satisfied. Indeed: 



+oo 



< /miijj^oo — / k(y, s, x, t)v(x, t)xK(x)dx < 



oo 



— [An(t- s)}- 1 ' 2 exp\M{t- s)}- (53) 

Cl 

[sup x€K v(x, t)) Hindoo exp(c 2 y 2 ) exp[- — - ^ dx] = 
if t — s > e for sufficiently small e > (like for example e = I/I6C2). 

It is our purpose to complete the previous analysis by demonstrating that, with 
the above assumption on v(x, t), the continuous Markov process we have constructed 
actually is the diffusion process. 



Our subsequent arguments will rely on the Dynkin treatise |17|. It is well known 
that the infinitesimal (local) characteristics of a continuous Markov process can 
be defined in terms of its, so called, characteristic operator. It is closely linked 
with the standard infinitesimal (Markov) generator of the process, and we shall 
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take advantage of this link in below. Let us agree, following Dynkin, to call a 
continuous Markov process a diffusion, if its characteristic operator U is defined on 
twice differentiable functions (we skip more detailed definition, JT7J). In this case 
and ) allow for the definition of a drift and diffusion function 

respectively: 

[U (x - x )) (x , s) = b(x , s) (54) 
[U((x - x ) 2 )}(x , s) = a(x , s) 

By results of Sections II. 3 and II. 4 we know that our transition probability density 
p(y,s,x,t) = k(y, s,x, t) 1^'% inspired by the Schrodinger boundary data problem, 
gives rise to a continuous Markov process. To see whether it can be regarded as a 
diffusion, we must verify the above two defining properties (54). 

At first, let us consider the infinitesimal operator A (Markov generator) of the 
corresponding strongly continuous semigroup T* : C OQ (R) — > C^R), which we have 
introduced via the formula (26). We are interested in domain properties of A, in 
view of the fact that the characteristic operator U is a natural extension of A, A C U, 

H • 

We denote C^(R) the space of continuous functions with compact support which 
possess continuous derivatives up to second order. For h G C^(R) we have 

1 r+oo 

lim &i0 -[ p(y,s,x,s + 5)h(x)dx - h(y)] = (55) 

J-oo 
1 1 f+°° 

— rlimsio -[ / k(y, s,x,s + S)v(x, s + 8)h(x)dx - v(y, s)h(y)\ 

v[y, s) o J-oo 

Because v is continuously differentiable with respect to time, we have 

v (x, s + 5) — v(x, s) + 5 d s v(x, s') (56) 

where (cf. the standard Taylor expansion formula) s' = s + fM, < i? < 1. Hence 

1 r +oc 

Hm sl0 -[ p(y,s,x,s + 6)h(x)-h(y))= (57) 

J-oo 
1 1 f+°° 

— rlims\ -[ / dx k(y, s,x,s + S)v(x, s)h(x) — v(y, s)h(x)\ + 

v{y, s) d J-oo 

\ r+oo 

— -limsi [ / Hy, s, x, s + 5)d s v(x, s')h(x)dx] 

v(y, s) J-oo 

We shall exploit the strongly continuous semigroup evolution associated with the 
parabolic system (5). Because of the domain property: C£°(R) C D(H) the smooth 



17 



functions with compact support are acted upon by H = A — c(x, s) and H is closed 
as an operator on C^R). But then also C%(R) C D(H) and so the first term in 
(57) takes the form: 

! r A(vh)(y,s)-c(y,s)v(y,s)h(y)] (58) 



v(y, s) 
while the second equals 



[d s v(y, s)]f(y) = -^—[-Av(y, s) + c(y, s)v(y, s)]f(y) (59) 



v(y,s) ' v(y,s) 
Thus, (55) is point-wise convergent: 

1 r+oo 



I r+oo 

limsiofi p(y,s,x,s + 5)h(x)dx - h(y)} = (60) 

J —oo 

[(Av(y, s)h(y) + 2Vv(y, s)Vh(y) + v(y, s)Ah(y) - c(y, s)v(y, s)h(y)- 



v(y,s) 



Vv 

(Av(y, s))h(y) + c(y, s)v(y, s)h(y)} = Ah(y) + 2( — )(y, s)Vh(y) 
Now, we shall establish the boundedness of: 

1 r+oo 

sup yeR . i0<s< e b I / P(y, s,x,s + S)h(x)dx - h(y) |] (61) 

J —co 

for some small e. 

Because C^(R) C D(H), so there holds 



\ r+oo 

-[/ k(y,s,x,s + 5)v(x,s)h(x)dx - v(y,s)h(y)} -> [A - c(y, s) (62) 

J-oo 



uniformly in as 5 — > 0. It implies that for any compact set K there is 

1 r+oo 

sup yeK ;o<s<e t\ p(y,s,x,s + 5)h(x)-h(y)\< (63) 

J —oo 

I I r+oo 

[sup y£K — A sup yeK . 0< s <e [-| / k(y, s, x, s + S)v(x, s)h(x)dx - v(y,s)h(y)\ + 

vyy, s) o J —oo 

/+oo 
k(y, s, x, s + 5)d s v(x, s')h(x)dx\] < oo 
-oo 

We have thus the required boundedness for all y E K i.e. on compact sets. 
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For y G R\K we shall make the following estimations. Because the support 
of h is compact, we can define supp h C [— n, n] for some natural number n. Let 
K = [Sn, 3n]. Then: 

sup yeR \ K;0<s< e -z | / p(y, s, x, s + 5)h(x)dx - h(y) | = 

J —oo 

sup y eR\K-,o<s<e^\ J p(y , s , x, s + 6)h(x)dx\ < (64) 
1 1 f 

s u P*6Jr ™p»efl\if!0<«<e 7-; r / k(y,s,x,s + S)v(x,s + S)dx < 

v(y, s) Jk 

1 1 r 

[sup xeK \h(x)\] [sup xGK . s < s >< s+e v(x,s')} sup y< z R \ K . fi< s <e ^ j k(y,s,x,s + 5)dx 

In view of our assumption v(y, s) > C\txp{— C2y 2 ), there holds: 

sup y eR\K-,o<s<e ^1 J p(y, s,x,s + 5)h(x)dx\ < (65) 

C ■ sup\ y \> 3n;0< s <£ exp(c 2 y 2 ) (T 3/2 J exp[ — — dx 

where 

C = Ci(47r)- 1/2 exp(Me) [sup xeK \h(x)\] sup xeK;s<s > <s+£ v(x, s') (66) 
If we choose e = I/I6C2, then 

exp(c2y 2 ) / exp[ — ] dx < 45exp( — — ) (67) 

J 4o 

for every \y\ > 3n, and so 

X f Ti^ 
sv>Py<zR\K;o<5<t ^ I J p{u, s, x, s + 5)h(x)dx\ < 4Csup 0<s<€ 5~ 1/2 exp( — — ) < 00 

(68) 

Consequently, the desired boundedness (62) holds true for all y G R, together with 
the previously established point-wise convergence (61). 

Altogether, it means, |]17|| , that the weak generator of T\ is defined at least on 
C 2 (R). Moreover, while acting on h G C 2 (R) it gives Ah + (Vlnv)Vh. Because T* 
is strongly continuous in C^R), the Markov generator A coincides with the weak 
generator, 0, i.e. A = A + (Vln v)V on C 2 (R). 

Finally, let us choose h G C 2 (R) such that h (x) = 1 in some neighbourhood 
of the point xq. Then, (x — xo)h (x) and (x — xo) 2 h (x) both belong to C 2 (R) and 
therefore: 

A[(x - x )h }(x , s) = A[(x - x )ho](x )+ (69) 
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2(Vlnv)(xo, s)V[(x — x )h Q ](xo) = 2(Vlnv)(x , s) 

A[(x - x ) 2 h }(x ,s) = 2 

Because A (Z U and U is a local operator, JlT), we have the following inclusion 
Cl(R) C D(U) and (we can get rid of h ): 

[U(x - x )](x , s) =2(Vlnv)(x ,s) (70) 

[U{x-x ) 2 ]{x ,s) = 2 

It means that we indeed obtain a diffusion process with the drift Vln v and a 
constant diffusion coefficient, according to the standards of || |27|, |2~8fl . 

It is worth emphasizing that since (x — xo)ho(x) and (x — xo) 2 ho(x) belong to 
D(A), and since functions from C 2 (R) can be used to approximate, under an integral, 
an indicator function of the set [xo — e, x$ + e], e > 0, we can directly evaluate: 

1 r+oa 

lim tis - / p(x ,s,x,t)(x - x )h (x)dx = (71) 



't-sJ- 

limti s / p(x , s,x,t)(x — x )dx = 2(Vln v)(x , s) 

t — S J \x—xo\<e 



and similarly 



lirritis / p(xo,s,x,t)(x — xq) dx = 2 (72) 

t — S J \x—xo\<e 

Because the Dynkin condition (28) implies that 

limtis / pfao, s, x, t)dx = (73) 

t S J \x — X()\>E 

we arrive at the commonly accepted definition of the diffusion process, summarized 
in formulas (71)-(73), with the functional expression for the drift, (71), given in the 
familiar, || ^7| f|, gradient form. 



III. Nonstationary Schrodinger dynamics: from 
the Feynman-Kac kernel to diffusion process 

In our previous paper 0, the major conclusion was that in order to give a definitive 
probabilistic description of the quantum dynamics as a unique diffusion process solv- 
ing Schrodinger's interpolation problem, a suitable Feynman-Kac semigroup must 
be singled out. Let us point out that the measure preserving dynamics, permitted 
in the presence of conservative force fields, was investigated in 0, see also [f2£], |T2|j . 
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The present analysis was performed quite generally and extends to the dynam- 
ics affected by time dependent external potentials, with no clear-cut discrimination 
between the nonequilibrium statistical physics and essentially quantum evolutions. 
The formalism of Section II encompasses both groups of problems. Presently, we 
shall restrict our discussion to the free Schrodinger picture quantum dynamics. Fol- 
lowing Ref. we shall discuss the rescaled problem so as to eliminate all dimen- 
sional constants. 

The free Schrodinger evolution idtip = —Aip implies the following propagation 
of a specific Gaussian wave packet: 

Mx,0) = (2*)- 1 f'exp(~) — (74) 



(c> i o-mA -1/2 .„f X 



2 



${x,t) = (-) 1/4 (2 + 2it)~ 1/2 exp[ 



7T 4(l + it) J 



So that 



p (x) = mx,0)\ 2 = (27T)- 1 / 2 exp[~] — > (75) 



p{x,t) = \4>{x,t)\ 2 = [2tt(1 +t 2 )]- 1/2 exp[- 



x 2 



2(l+t 2 ) J 

and the Fokker-Planck equation (easily derivable from the standard continuity equa- 
tion dtp = — V(vp), v(x,t) = xt/(l + t 2 )) holds true: 

d t p = Ap- V(6p) , 6(0;, t) = x (76) 

1 + V 

The Madelung factorization ip = exp(R + iS) implies (notice that v = 2V5* and 
b = 2V(i? + S)) that the related real functions 9(x,t) = exp[R(x,t) + S(x,t)] and 
9*(x,t) = exp[R(x,t) — S(x,t)] read: 

x 2 1 — t 1 

6(x,t) = [2tt(1 + t 2 )]' 1/i exp(-— Y^~p ~ 2 arctant " > 

0*{x,t) = [27r(l + t 2 )]-V4 exp (_^ 1±1 + l rc£anf ) ( 77 ) 

A 1 + t z 2 

They solve a suitable version of the general parabolic equations (5), namely : 

d t 9 = -A9 + c6 (78) 
d t 0* = AO* - c9* 
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with 



,2 i A 1/2 



t(l - t,= 2(rw-iTf =2 ^ (79) 

Anticipating further discussion, let us mention that the Feynman-Kac kernel, 
in this case, is a fundamental solution of the time adjoint system (78). For clar- 
ity of exposition, let us recall that a fundamental solution of the forward parabolic 
equation (5) is a continuous function k(y,s,x,t), defined for all x, y, G R and all 
< s < t < T, which has the following two properties: 

(a) for any fixed (y, s) e Rx (0, T), the function (x, t) — > k(y, s, x, t) is a regular (i.e. 
continuous and continuously differentiable the needed number of times) solution of 
the forward equation (5) in R x (s, T] 

(b) for all continuous functions 4>(x) with a compact support, there holds lim(t tX )^(s,z) 
/-« s > t)(j)(y)dy = (f)(z). 

First, we need to verify (this will be done self-explanatorily) that c(x,t), (79), is 
Holder continuous of exponent one on every compact subset of R x [0, T]. It follows 
from direct estimates: 

\^M)-c{x lM )\<\\j^ 2 -j^ 2 \ + li^-i^l^ (80) 
!, x 2 Xi \x 2 \ , , , |*i| + |t 2 | 

+ t2 — tl 



2'l+t| l+tf'V 1 + t 2 1+t 2/ I- il (l + f 2)( 1 + f2) 

But, in case of |xi|, \x 2 \ < K and \t\\, \t 2 \ <T we have 



\c{x 2 M)-<x 1 M)\<K\-^--^\ + 2T|i 2 -* 1 | (81) 

1 + t 2 1 ~T t± 



Furthermore: 



I x 2 x l i ^ i ^2 — x l i i ^2^1 — ^1^2 I ^ /oo\ 

1 T+1 "l + t! 1 " 1 (i + ti)(i + t?)' + '(l + tixi + t?) 1 " ( } 

|x 2 -x 1 \+ T 2 \x 2 - Xl \ + 2irT|t 2 - tx| 
implies (the new constant C majorizes all remaining ones) 

\c(x 2 ,t 2 ) -c(xi,fi)| < C(|or 2 + |*2 — *x|) < v^C [(x 2 -xi) 2 + (t 2 -ti)T /2 

(83) 

Let us also notice that we can introduce an auxiliary function h(x, t) = arctan t 
such that there holds 

Ah - c(x, t)h - d t h = ~ 0^r 2 < (84) 
22 



We have thus satisfied the crucial assumptions I and II of Ref. |2~5|| . As a conse- 
quence, we have granted the existence of a fundamental solution k(y,s,x,t) > 0. 
Moreover, for every bounded and continuous function 0(x), |0(x)| < C, where C > 
is arbitrary, the function 

/+oo 
k(y,0,x,t)(f>(y)dy (85) 
-oo 

is a solution of the Cauchy problem, i.e. solves (79) under the initial condition 
u(x,0) = (f)(x), so that \u(x, t)\ < C. All that implies the uniqueness of the fun- 
damental solution k(y,s,x,t), and in view of —c(x,t) < 1 its strict positivity. The 
function k(y,s,x,t) is also a solution of the adjoint equation with respect to vari- 
ables y, s: d s k = —A y k + c(y,s)k in R x [0,T). It is obvious that the Chapman- 
Kolmogorov composition rule holds true, in view of the validity of the Feynman-Kac 
representation in the present case. 

Basically, we must be satisfied with the Feynman-Kac representation of the fun- 
damental solution, whose existence we have granted so far. In our case, the so 
called parametrix method, ||21|| , can be used to construct fundamental solutions. In 
fact, since c(x,t) is locally Lipschitz i.e. Holder continuous of exponent one and 
quadratically bounded \c(x, t)\ < x 2 + 1, the infinite series: 

+oo 

k(y,s,x,t)= £(-l) n A; n (y, S ,:r,t) (86) 

n=0 

where k (y, s, x, t) = [An(t — s)}^ 1 / 2 exp[— (x — y) 2 /4(£ — s)] is the heat kernel and 

rt r+oo 

k n (y,s,x,t)= dr dz c(z,r) k n -i(y, s, z,r)k (z,T,x,t) (87) 

J S J —CO 

are known to converge for all x, y G R, < s < t < T, and t — s < T where T < T, 
and define the fundamental solution, |30| . 

By putting p(y,s,x,t) = k(y, s, x, t) ffi '*( we arrive at the fundamental solution 
of the second Kolmogorov (Fokker-Planck) equation 

dtp(y, s, x, t) = A x p(y, s, x, t) - V x [b(x, t)p(y, s, x, t)) (88) 

where 6 = 2^ and p = uv, and in particular p = 99* = \ip\ 2 , are consistently 
propagated by p. It is the transition probability density of the Nelson diffusion 
associated with the solution (74) of the Schrodinger equation, and at the same time 
a solution of the first Kolmogorov (backward diffusion) equation 

d s p(y, s, x, t) = -A y p(y, s, x, t) - b(y, s)W y p(y } s, x, t) (89) 
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Equations (88), (89) prove that the pertinent process is a diffusion: it has the 
standard local (infinitesimal) characteristics of the diffusion process, |2(J . 

Obviously, the above definition of p in terms of k induces the validity of the 
compatibility condition 

c(x,t) = 2[d t ln9(x,t) + ^[^^ + Vb(x,t)} (90) 

connecting the drift of the diffusion process with the Feynman-Kac potential gov- 
erning its local dynamics: cf. Refs. p9| , |J and |31|] where the Ehrenfest theorem 
analogue was formulated for general (non-quantal included) Markovian diffusions. 

Let us point out that our quantally motivated example was chosen not to show 
up a typical for quantum wave functions property of vanishing somewhere. In fact, 
because of restricting our considerations to strictly positive Feynman-Kac kernels 
and emphasizing the uniqueness of solutions, we have left aside an important group 
of topics pertaining to solution of the Schrodinger boundary data problem when: 

(i) the boundary densities have zeros 

(ii) the interpolation itself is capable of producing zeros of the probability density, 
even if the boundary ones have none. 

Only the case (i) can be (locally) addressed by means of strictly positive semigroup 
kernels, however the uniqueness of solution is generally lost in space dimension higher 
than one, 0, |[ General existence theorems are available p8 , p9| and indicate 



that one deals with diffusion-type processes in this case, see e.g also 0, [§, 



The case (ii) seems not to be ever considered in the literature, see however 32 
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